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The mass profile of A1413 observed with XMM-Newton: 
implications for the M-T relation. 
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Abstract. We present an XMM-Newton observation of A1413, a hot (kT = 6.5 keV) galaxy cluster at z = 0.143. 
We construct gas and temperature profiles over the radial range up to ~ 1700 kpc. This radius corresponds 
to a density contrast 5 ~ 500 with respect to the critical density at the redshift of the cluster, or equivalently 
~ 0.7r2oo- The gas distribution is well described by a (3 model in the outer regions, but is more concentrated in the 
inner ~ 250 kpc. We introduce a new parameterisation for the inner regions, which allows a steeper gas density 
distribution. The radial temperature profile does not exhibit a sharp drop, but rather declines gradually towards 
the outer regions, by ~ 20% between 0.1r2oo and 0.5r2oo- The projected temperature profile is well described by 
a polytropic model with 7 = 1.07 ±0.01. We find that neither projection nor PSF effects change substantially the 
form of the temperature profile. Assuming hydrostatic equilibrium and spherical symmetry, we use the observed 
temperature profile and the new parametric form for the gas density profile to p roduc e the total mass distribution 
of the cluster. The mass profile is remarkably well fitted with the Moore et al. (1999) parameterisation, implying 
a very centrally peaked matter distribution. The concentration parameter is in the range expected from numerical 
simulations. There are several indications that beyond a density contrast 5 ~ 600, the gas may no longer be 
in hydrostatic equilibrium. There is an offset with respect to adiabatic numerical simulations in the virialised 
part of the cluster, in the sense that the predicted mass for the cluster temperature is ~ 40% too high. The gas 
distribution is peaked in the centre primarily as a result of the cusp in the dark matter profile. The X-ray gas to 
total mass ratio rises with increasing radius to / gas ~ 0.2. These data strongly support the validity of the current 
approach for the modeling of the dark matter collapse, but confirm that understanding the gas specific physics is 
essential. 

Key words. Galaxies: clusters: individual: A1413, Galaxies: clusters: Intergalactic medium, Cosmology: observa- 
tions, Cosmology: dark matter, X-rays: galaxies: clusters 



1. Introduction 

The simple first-order formation scenario for galaxy clus- 
ters, in which they grow through the gravitational infall 
and subsequent merging of smaller subunits, provides a 
remarkably good description of the large-scale properties 
of these objects. Within this hierarchical model, the gas 
trapped in the potential well of a cluster is heated to the 
observed X-ray emitting temperatures by the shocks due 
to the formation process; merger features in the gas dis- 
tribution are then erased in roughly a sound crossing time 
(~ few Gyr), leaving the gas in hydrostatic equilibrium 
(HE). 

Observation of this gas is a powerful tool for un- 
covering the physical characteristics and formation his- 
tory of a cluster. Substructure in X-ray images, combined 
with optical data, can give clues to the dynamical state 
(e.g., Buote 2001). Direct (temperature maps) and indi- 



rect (hardness ratio maps) methods can give an indication 
of where (and if) interactions and mergers are still occur- 
ring (e.g., Markevitch et al. 1999; Neumann et al. 2001). 



In addition, for clusters in reasonably relaxed state, the as- 
sumption of HE and spherical symmetry allow the deriva- 
tion of the spatial distribution of both the gas and total 
cluster mass by using the information from the X-ray sur- 
face brightness and temperature profiles. This approach, 
which is of fundamental use in cluster studies, has been 
shown to give masses which are accurate to about ±20% 
when applied to simulated clusters (e.g., Evrard, Metzler 
& Navarro |l996| [EMN96] ; Schindler pL996| ). 
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Numerical simulations based on gravitational collapse 
are an essential counterpoint to the observations, being as 
they are ideal scenarios with exactly measurable quan- 
tities, thus offering a direct comparison with the real 
data. A crucial result from these simulations is the sug- 
gestion that CDM haloes with masses spanning several 
orders of magnitude follow a universal density profile in- 
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dependent of halo mass or cosmology (Navarro, Frenk & 
White |1997| [NFW]). As the X-ray emitting gas lies in 
the potential well of the CDM halo, this suggests that 
many directly measurable cluster properties should dis- 
play self-similarity. This is observationally testable and 
indeed, regularity in the local cluster population has been 
found in previous ROSAT, ASCA and BeppoSAX stud- 
ies, where the gas density and temperature profiles of hot, 
relaxed clusters do appear similar when scaled to units 
of the virial radius^ (Markevitch et al. 1998[ Neu mann 
& Arnaud [199$ V ikhlinin, For man & Jones |1999|: Irwin 
& Bregman 120001; De Grandi & Molendi |2002[ Arnaud, 



Aghanim & Neumann 2002). The very existence of these 
similarities gives strong support to an underlying univer- 
sality in the dark matter distribution, leading to a pleasing 
convergence between the observed and simulated proper- 
ties of galaxy clusters. 

However, the temperature profiles in particular have 
generated much discussion, as rather different profile 
shapes have been found for similar samples observed 

by 



1998 



the same satellite (e.g., Markevitch et al 
White p000|[^ gC^]; Irwin & Bregman pOOOj De Grandi & 
Molendi 2002 [BeppoSAX]). These studies have been ham- 
pered somewhat by both PSF issues and sensitivity limits. 
The former has an inevitable effect on the spatial resolu- 
tion and is a possible source of systematic uncertainty, the 
derivation of the profiles being potentially sensitive to the 
exact correction for the PSF and the detailed modelling 
of the non-resolved cooling flow component. The latter 
leads to an inability really to constrain parameters be- 
yond the supposedly isothermal regime, which is expected, 
from simulations, to extend to ~ 0.5r2oo- As a direct con- 
sequence of this, there are relatively few galaxy clusters 
for which sufficiently high quality data were available for 
an accurate determination of the total mass and the cor- 
responding density profile. Furthermore, any systematic 
uncertainty in the shape of the radial temperature distri- 
bution can have a direct effect on the derived mass. For 
example, the temperature profile obtained by Markevitch 



et al. (1998) gives mass values that are 1.35 and 0.7 times 
that derived assuming isothermality at 1 and 6 core radii 
respectively. As a result, the actual form of the density 
profile is still a largely untested quantity, at least from an 
observational point of view. 

Clusters can also be used to provide cosmological con- 
straints. For any given cosmology and initial density fluc- 
tuation, the mass distribution of virialized objects can be 
predicted for any given redshift. Constraints on cosmolog- 
ical parameters, as and O, can be found by comparing the 
predictions with the observed cluster mass function and its 
evolution (Perrenod 198C). For this, however, a great num- 
ber of accurate observational masses are needed. In the 
calculation of the observed cluster mass function, the stan- 
dard way to overcome the paucity of data is to use average 



cluster temperatures, taking advantage of the tight mass- 
temperature relation predicted by numerical simulations, 
where M oc T 3 / 2 (e.g., EMN96). While observations have, 
for hot clusters at least, recovered the slope of this rela- 
tion, observed masses imply a normalisation consistently 
lower than found by simulations (e.g., Horner, Mushotzky 
& Scharf |1999|; Nevalainen, Markevitch & Forman |2000: 



1 Normally defined from numerical simulations as the radius 
of fixed de nsity c ontrast S — 200, or r2oo (e.g., Evrard, Metzler 
& Navarro |l996fl 



Finoguenov, Reiprich & Bohringer 2001). However, these 
total cluster mass estimates, except in a few cases, re- 
quired an extrapolation of the data and the level of the dis- 
crepancy is sensitive to the assumed temperature profile 
(e.g. see Horner et al. 1999, Neumann & Arnaud 1999). 

XMM-Newton and Chandra offer, for the first time, 
sufficiently good spatial and spectral resolution for self- 
consistent determinations of global cluster observables 
such as gas density, temperature and mass profiles. We 
are now observing clusters with unsurpassed clarity. 
Chandra, with higher resolution, is the instrument best- 
suited for the study of cluster cores. In the most recent 
Chandra study by Allen, Schmidt & Fabian (2001a), mass- 
temperature data from 6 clusters are measured up to r25oo, 
and compared to the reference simulations of EMN96 and 
Mathiesen & Evrard ( |2001| ). Once again, a systematic off- 
set of ~ 40% is found between the observed and simulated 
M — T curves, in the sense that the predicted tempera- 
tures are too low for a given mass. XMM-Newton, with 
its high throughput and large field of view, is the satellite 
best-matched for the study of the larger scale structure of 
these objects, and for the determination of essential quan- 
tities out to a good fraction of the virial radius. With this 
capability it is possible to test for other effects, such as 
potential variations of the normalisation with radius. 

In this paper, we use XMM-Newton observations of the 
relaxed cluster A1413 at z = 0.143 to derive the large scale 
properties to high resolution, and compare the results to 
those obtained from both observations and simulations. 
We address several questions which have been the subject 
of a large amount of debate in the literature. In particular, 
we compare our temperature profile with previously de- 
rived composite profiles from large samples observed with 
ASCA and BeppoSAX , and we compare both the form 
and normalisation of our mass profile with that expected 
from numerical simulations. 

Throughout this paper we use Ho = 50 km s _1 Mpc -1 , 
and unless otherwise stated, ft m = 1 and fl\ — (go — 
0.5). In this cosmology, at the cluster redshift of z = 0.143, 
one arcminute corresponds to 196 kpc. 



2. Data analysis and preparation 

2.1. Observations 

A1413 was observed in Guaranteed Time for 29.4 ks dur- 
ing XMM-Newton revolution 182 (2000 December 16). 
Calibrated event files were provided by the XMM-Newton 
SOC. The MOS and pn data were obtained with the THIN1 
and MEDIUM filters, respectively. For the pn data set, we 
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extracted single events, which correspond to PATTERN 0, 
while for the MOS data sets PATTERNS 0-12 were selected. 

Dedicated blank-sky data sets, which consist of sev- 
eral high-galactic latitude pointings with sources removed 
(Lumb 2002 ), were used as background for the whole 
of this analysis. These data sets are distributed as cal- 
ibrated event files which have already been treated with 
the SAS. We extract the background events using the same 
PATTERN selection criteria as outlined above. In addition, 
we transformed the coordinates of the background file such 
that they were the same as for the A1413 data set. In this 
way we can ensure that all source/background products 
come from exactly the same regions of the detector, thus 
minimising detector variations. 



2.2. Vignetting correction 



The method described in Arnaud et al. (2001b) was used 
to correct spectra and surface brightness profiles for vi- 
gnetting effects. Briefly, this method involves weighting 
each photon with energy (E), detected at position (xj, i/j), 
by the ratio of the effective area at the detected position 
A Xj , yj {E) to the central effective area A 0y0 (E). 

The background data were treated in the same man- 
ner as the source. Note that the background component 
induced by cosmic rays (see below) is not vignetted, but 
since source and background observations are treated in 
the same way, the correction factor is the same and thus 
cancels. 



2.3. Background subtraction 

The XMM-Newton background, consisting of several com- 
ponents, is both time and energy-dependent, and so sub- 
traction is a subtle process. Furthermore it is essential that 
the subtraction is done correctly, especially so for extended 
sources like clusters of galaxies, where background effects 
begin to play a role at large off-centre distances where 
the surface brightness declines approximately as r -4 (e.g., 



Vikhlinin et al. |1999| ). 

The soft proton background due to solar flares can- 
not be corrected for in the normal fashion (e.g., spectral 
subtraction) as it displays extreme temporal and flux vari- 
ability, causing the spectrum to change rapidly with time. 
At the moment it can only be removed by excising all 
frames above a certain count-rate threshold, the main ef- 
fect of which is to considerably reduce the effective ex- 
posure time. For these observations, the 3er threshold for 
each camera was calculated using the method described in 
Appendix A, and all frames not meeting this criterion were 
rejected. In practice the observation is very clean. Note 
however that the pn is considerably more sensitive to the 
flares. The final exposure times were 24163s, 24567s and 
10254s for MOS1, MOS2 and pn cameras, respectively. 
The blank-sky backgrounds were cleaned using the same 
criteria. 



The blank-sky background represents effectively the 
particle induced background, dominant in the hard X-ray 
band, which is, both spatially and temporally, relatively 
constant. Nevertheless, this background is variable at the 
~ 10%-level, and so it is frequently necessary to normalise 
the background. We normalise these observations using 
the count rate in the [10-12] keV and [12-14] keV bands, 
for MOS and pn respectively, treating each camera sep- 
arately. We varied the normalisation by ±10% to assess 
any systematic uncertainties. 

However, the blank-sky data set does not necessar- 
ily represent the cosmic X-ray background (CXB), be- 
cause this is variable across the sky, espec ially the soft 
X-ray component (see Snowden et al. 1997 ). We use th e 
method described in P ratt, Arnaud & Aghanim ( 2001 ), 
and Arnaud et al. ( 2002] ) to correct for the difference of the 
CXB. An annular region external to the cluster emission 
(between 9' and 13' in this case) is used to estimate the 
local background. The normalised spectrum of the same 
region of the blank-sky background is then subtracted, 
giving a difference spectrum, which can then be scaled ac- 
cording to the size of any extraction region and subtracted 
directly from the source spectra. A similar procedure is 
applied to subtract the residual CXB component for the 
surface brightness profile (see Arnaud et al. 2002 for de- 
tails). 

In addition to the above, the pn data were corrected for 
'out of time events', which occur when a photon hits the 
CCD during the read-out process in the imaging mode. 

3. Morphology 

3.1. Image 

We show in Figure^ the MOS1+MOS2 image of the cluster, 
produced simply by adding the data from each camera 
without accounting for vignetting. The image is striking: 
the cluster displays an unmistakably elliptical shape, with 
a clear brightness enhancement directly to the south, and 
there are a large number of sources in the field of view. 

3.2. 2D (3-model fitting 

Motivated by the apparent excess of counts to the south 
of the cluster (see Figure [TJ) we fitted the MOS1 +MOS2 
image with a 2D /3-model in order to quantify the signif- 
icance of this feature. In fitting the image, we followed 
closely the procedure described in Neumann & Bohringer 



(1997). Images were extracted in the [0.3-1.4] keV band 



from the weighted MOS event files in pixels of size 3". 4 
and added to make a combined MOS image. Since in the 
case of weighted events Poissonian errors do not apply, 
errors were correctly calculated from the weights using 



• 2 (see Arnaud et al. 2001b). An error image 

was generated for each instrument, and these images were 
added quadratically. The fitting procedure described be- 
low was tested and optimised on simulated data before 
application to the real data. 



XMM-Newton observation of A1413 




Table 1. 2D /3-model fits, 1'.3 - 13', MOS image 



11" 56 



Rfqht Ascension 



Fig. 1. The MOS1 +MOS2 counts image of the whole field 
of view of the A1413 observation. Note the large number of 
sources and in particular, what appears to the an extended 
source to the south of the cluster itself. 



The x 2 test used in the fitting procedure assumes 
Gaussian statistics, for which the mean is the most likely 
value. The statistics are actually not Gaussian in the ex- 
ternal regions of the field of view, dominated by the back- 
ground. In these regions the number of photons per pixel 
is low and follows a Poisson distribution for which the 
mean is larger than the most likely value. If the image is 
not smoothed before fitting, there is thus a tendency to 
underestimate the mean background level, leading to er- 
roneous values for the fitted cluster parameters. The com- 
bined MOS image was thus smoothed with a Gaussian of 
with a = 10" before fitting. The error image was treated 
according to the error propagation function for Gaussian 
filtered images described in Neumann & Bohringer ( 1997 ). 
We fix all error pixels with a value of to have a value 
of 1 before fitting, meaning that we can use \ 2 fitting but 
are unable to determine confidence parameters on the fit. 
The data were then fitted with a 2D /3-model of the form: 



S(x,y) = S Q (l + F 1 +F 2 ) 



-3/3H 



B 



(1) 



where 



F, = 



[cosev(x — xq) + sina(y — yo)] 2 
[— sina(x — Xo) + cosev(y — yo)Y 



Here, xo, yo is the position of the centre of the cluster; x, y 
are the coordinate positions of each pixel; a\ , a 2 are the 
major and minor core radii; a is the position angle; and 
the background is included in the model via the parameter 
B. 



Parameter 



P 

r c long 
r c short 
PA 

Centre a 
Centre Dec 



0.72 

284.6 kpc 
201.2 kpc 
2°26' 

ll h 55 m 18 s .9 
23°24'13".8 




Fig. 2. Residuals after subtraction of the 2D /3-model, 
smoothed with a Gaussian with a — 5/V2 pixels (~ 12"). 
This is a zoomed image where the centre of the cluster is the 
bright elliptical region at the centre of the image (a possible 
cooling flow?). The dynamic range is from -10 to +20cr, where 
areas of low a are black and areas of high a are white. Contours 
are between +2 and +10ct, in steps of 2. 



We fit the image between 1'.3 (see below) and 13' from 
the cluster center, excluding obvious point sources. The re- 
sults of the 2D fit are shown in Table [j]. Note that the fit- 
ted parameters are slightly dependent on the outer radius 
and the a of the Gauss filter, but the results are always in 
good agreement with the ID fit, described below. 

In order to quantify the significance of the excess of 
flux to the south of the cluster, we subtract the 2D j3- 
model from the data and calculate the significance of the 
residuals using the prescripti on de scribed in the Appendix 
of Neumann & Bohringer (|1997|) . The excess is an ex- 
tended source detected at > lOer, as shown in Figure 0. 

We extracted the spectrum from a circular region of 
radius ~ 1' centred on the excess. This spectrum unfor- 
tunately does not contain sufficiently strong line emission 
for a redshift estimate, so we fitted using a MEKAL model 
with the same redshift as A1413, absorbed with the galac- 
tic column density toward the cluster (2.19 x 10 20 cm~ 2 
from Dickey & Lockman |1990p . We find kT = 3.1 keV. 
An overlay of the significance contours on the DSS plate 
of the image did not reveal any obvious sources associ- 
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ated with the excess, and a hardness ratio map did not 
reveal any interaction with the main cluster. Our tenta- 
tive conclusion is that the source is either a foreground 
or background cluster: deeper optical observations of the 
region should resolve the issue. 

4. Gas density profile 

4.1. Surface brightness profile 

For each camera, we generated an azimuthally averaged 
surface brightness profile for both source and background 
observations. Weighted events from the corresponding 
event files were binned into circular annuli centred on the 
position of the cluster emission peak. We cut out serendip- 
itous sources in the field of view and the southern sub- 
structure. The background subtraction was performed as 
described in Sec. 0. We consider the profiles in several en- 
ergy bands. Due to the contribution of the instrumental Al 
K line around 1.5 keV, we ignored the [1.4-2.0] keV band 
in both cameras. To maximise the signal to noise (S/N) 
ratio, particularly in the outer cluster region, we choose to 
base the following on analysis of the [0.3-1.4] keV band. 

We checked that the vignetting corrected and back- 
ground subtracted profiles of the three cameras are con- 
sistent: they differ only by a normalisation factor within 
the error bars. We thus coadd the profiles and bin the 
resulting profile is the following way. Starting from the 
central annulus, we re-binned the data in adjacent annuli 
so that i) at least a S/N ratio of 3a is reached and ii) the 
width of the bin increases with radius, with A(0) > 0.18. 
Such a logarithmic radial binning insures a S/N ratio in 
each bin roughly constant in the outer part of the profile, 
when the background can still be neglected. 

The resulting surface brightness profile, S(8), is shown 
in Fig. |[ The cluster emission is significantly detected up 
to R det = 8.6' or 1.7 Mpc. 

4.2. Density profile modelling 

We fitted S(9) with various parametric mode ls convolved 
with th e XM M-Newton PSF (Ghizzardi pOOl] Griffiths & 
Saxton 2002), binned into the same bins as the observed 
profile. 

A single /3-model cannot account for the data. When 
the entire radial range is fitted, the reduced \ 2 is ~ 13; for 
the best fit slope, (3 — 0.60, and core radius, 9 C — 0.66'. 
An excess of emission is readily apparent in the centre 
and a lower reduced \ 2 is obtained when excluding the 
central region from the fit. The reduced % 2 decreases with 
increasing cut-out radius until it stabilises for i? cu t ~ 1.3'. 
In that case we obtained \ 2 — 47 for 31 d.o.f, with j3 = 
0.71 ± 0.02 and 6 C = 1.30' ± 0.09'. The best fit model 
is plotted as a dotted line in Fig. ^. The /3 value is not 
surprisingly larger than the value (/3 = 0.62) derived by 
Cirimele, Nesci & Trevese (1997) from their global fit to 
the ROSAT profile, but is in excellent agreement with the 
value j3 = 0.70 ± 0.02 obtained by Vikhlinin, Forman & 
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Fig. 3. Combined MOS1, MOS2 and pn surface bright- 
ness profile of A1413 in the [0.3 — 1.4] keV energy band. 
The profile is background subtracted and corrected for vi- 
gnetting effects. Black (red) [green] lines: best fit KBB 
(Eq.g), BB (Eq.|) and AB (Eqg) models convolved with 
the XMM-Newton PSF binned as the observed profile. 
Dotted line: best fit /3-model fitted to the outer region of 
the cluster (9 > 1.3'). See Sec.[h^ for model details and 
Table || for best fit parameter values. 



Jones ( 1999 ) by fitting the outer cluster region {9 > 4.9'). 
There is also an excellent agreement between the ID and 
2D (3- values. 

We note that the last two points (7.75' < 9 < 8.6') lie 
significantly below the best fit model (a 3a effect for the 
last bin). The cluster flux in the last bin is about 16% of 
the total background and we cannot totally exclude that 
this discrepancy is an artifact due to remaining system- 
atic uncertainties in the background subtraction. This is 
further discussed in Sect. 9.4. These last two points are 



discarded in the present analysis. 

For the mass analysis which follows (Sect. ||) it is 
convenient to have an analytical description of the gas 
density radial profile (nn(r)) at all radii. We thus tried 
several alternative parameterisations, with behaviour at 
large radii similar to a /3-model: 

— AB model: A cusped profile similar to the NFW uni- 
versal density profile: 



n H (r) =A[ — 
r, 



1+ I - 

TV 



3ii 



(2) 
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where a is the slope at small radii. 



Table 2. Results of the surface brightness profile fits. 



— BB model: We also introduce a double isothermal 
model (BB), assuming that both the inner and outer 
gas density profile can be described by a /3-model, but 
with different parameters. 



Parameter 



AB model 



BB model 



KBB model 



r < R cut n H (r) = n H , 
r > R cut n H (r) = N 



3ft, 



(3) 



The boundary between the two regions, i? C ut, is a free 
parameter of the model and we took care that the den- 
sity distribution is continious across R C ut, as well as its 
gradient (for continuity of the total mass profile, see 
Eq. |l|): 

-3/3 in 

: L ^ Sr (4) 



N = n H .o 



and 



An = 



Rout 



Rout 



1 



(5) 



— KBB model: We finally consider a generalisation of 
the /3-model for the inner region, allowing a more cen- 
trally peaked gas density profile in the core: 



r < R cut n H (r) 

r > R cut n H (r) 
where £ < 1 and 

1 _l_ I E™ 



«h,o 



N 



1 



25 



, _M 



(6) 



N 



HH,0 



2 S C 



-3ft, 



1 



flout 



f 



An = 



2? 



(7) 



(8) 



The parameters £ and r Cj i n are strongly correlated. An 
arbitrary low value of £ can fit the data if no upper 
limit is put on r Cj j n . The lower limit on this parameter 
given in Table ^ is obtained by imposing r C: j n < 1'. 

The corresponding surface brightness profile is com- 
puted numerically by integration of the emission measure 
along the line of sight: 

n|(6) 



Sx(0) (x 



-.dr z 



(9) 



n H ,o (Kr'cnT^) 




2.15 


3.07 




1.54' 


1.29' ±0.10' 


1.34' ±0.12' 


/? 


0.69 


0.71 ±0.02 


0.71 ±0.02 






1.47' ±0.13' 


i.69'±8:l| 


fc,in 






0.41'" .i3 


Z 






0.49l°;^ a 


a 


0.68 






xVdof 

Xrcd 


112/51 


70.4/48 


64.8/47 


2.20 


1.47 


1.38 



Notes: All errors are at the 90% confidence level. 
a The maximum value of r c t n is fixed to 1'. 



where r = d^0 and is the angular distance. The emis- 
sivity in the considered energy band was estimated using 
an absorbed isothermal model at the cluster mean temper- 



ature (given Sect. 5.1), taking into account the instrument 
response. In the soft energy band considered, this emis- 
sivity is insensitive to the observed temperature gradient 
(shown Sect. |5.2| ). Note that the profile beyond R cu t ob- 
tained for the BB and KBB models is a classical /J-model. 
The inner surface brightness profile for the BB model can 
be analytically computed using incomplete Beta functions. 

The best fit models are plotted in Fig. ||, together with 
the residuals. The corresponding best fit parameters with 
errors and x 2 values are given in Table |^. In all cases, the 
outer slope, 0, is consistently found to be similar to the 
slope obtained by fitting only the outer part of the profile. 
We found that the AB model does not provide a partic- 
ularly good representation of the data: the reduced % 2 is 
Xrcd ~ ^ and the residual profile below i? cu t clearly indi- 
cates that the gas distribution is less peaked than a cusped 
profile. In other words, the gas distribution possesses a 
core. The best fit is obtained with the KBB model, but 
the reduced % 2 , x 2 c d ~ 1-38 is still larger than 1. However, 
the residuals are small (at the 3% level on average) and 
might be due in part to the observed departure from spher- 
ical symmetry. As the BB model is a special case of the 
KBB model (£ fixed to £ = 1), we can compare both 
models using a F-test. The KBB model provides a better 
fit than the BB model at the 95% confidence level, sug- 
gesting that the density distribution in the core is indeed 
more centrally peaked than for a conventional A -11101 ^!. 
This KBB model is thus adopted for the remainder of the 
analysis. 

5. Spatially resolved spectroscopy 

5.1. Global spectrum 

For each instrument, a global spectrum was extracted 
from all events lying within 5'.1 from the cluster emission 
peak. This radial range was chosen to maximise the 
S/N ratio, allowing us to check in detail the consistency 
between the three cameras. Each global spectrum was 
fitted with an absorbed MEKAL model with the red- 
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shift fixed at z — 0.143. The normalisation for each 
instrument was left as an additional free parameter. 
We excluded the energy bins around the strong fluo- 
rescence lines of Ni,Cu & Zn from the pn spectrum fit. 
These lines, present in the background, are not well 
subtracted by the procedure described in Sec. || because 
they do not scale perfectly with the continuum of the 
particle-induced background. In all fits we used the fol- 
lowing response matrices: ml_thinlv9q20t5r6_all_15.rsp 
(MOS1), m2.thinlv9q20t5r6_all.15.rsp (MOS2) and 
epn_ef20_sY9 .medium. rsp (pn). 

Fitting the data from all instruments above 0.3 keV, 
with the absorption fixed at the galactic value of TVh = 
we found inconsistent values for the 



Table 3. Influence of the low-energy cutoff. Absorption values 
in bold are frozen at the galactic value. 



2.2 x 10 20 cm 



temperature derived with the MOS and pn cameras: 
kT = 6.91±g;|| keV (MOS1), 6.33t°;^ keV (MOS2) and 
5.76lg;}g keV keV (pn). A better agreement between the 
cameras, together with a lower \ 2 value, are obtained if 
the TVh value is let free, but then the best fit TVh values are 
significantly lower than the 21 cm value (see Table ||). In 
other words the data presents an excess at low energy as 
compared to an isothermal model absorbed with the galac- 
tic hydrogen column density. This effect could be due to a 
true soft excess component (e.g Durret et al. 2002) and/or 
an artifact due to remaining calibration uncertainties. In 
particular, it is known that the EPIC-pn and MOS cam- 
eras show a relative flux difference which increases with 
energy above 4.5 keV, resulting in a MOS spectral slope 
flatter than the pn (Saxton 2002; Griffiths et al 



2002) 



We then performed a systematic study of the effect of 
imposing various high and low-energy cutoffs for each in- 
strument. The A^h is fixed to the 21 cm value. Having first 
found that progressive cutting of the high energy chan- 
nels had a negligible effect on the derived temperatures, 
we then varied the low energy cutoff, for which the results 
are shown in Table ||. 

This table shows that there is an optimum low-energy 
cutoff for each instrument, above which no amount of fur- 
ther cutting of the low-energy response will significantly 
affect the temperature. The temperature stabilises above 
a certain cutoff point for each instrument, this being ~ 0.6 
keV for the MOS cameras and ~ 1.0 keV for the pn 
camera. The adoption of these low-energy cutoffs has the 
pleasing effect of bringing the temperatures for each in- 
strument into agreement both with each other and with 
previous ASCA analysis. The combined MOS -fpn global 
temperature is kT = 6.85±g;^ keV (90% confidence for 



fi 85+ - 15 

°- oo -0.16 

one interesting parameter, \ 2 — 1459.6 for 1436 dof) in 



agreement with the results of Ikebe et al. ( 2002), wh o find 
kT = 6.56±g |5 keV and Matsumoto et al. ( |200l| ), who 
find kT = 6.72 ±0.26 keV, and marginally consistent with 
the result of White ( [2000D , who finds kT = 7.32±g-j| keV. 

It thus appears that the discrepancies observed by fit- 
ting the whole energy range are mostly due to some resid- 
ual calibration uncertainties in the low-energy response of 
all instruments and/or a true soft excess. The scientific 
analysis of such a possible soft excess is beyond the scope 



Instrument 


Band 


A H 


kT 


X7dcf 




(keV) 


(xlO 20 cm" 2 ) 


(keV) 




MOS1 


> 0.3 


1 04 +UM 


7 ci +0.40 
' ■ J1 — 30 


394.3/395 




> 0.3 


2.19 


fi qi+0-23 

v-i3±_Q 23 


424.2/396 




> 0.6 


2.19 


7 i r+0.25 
' ■ 1 ' J -0.25 


386.7/376 




> 0.8 


2.19 


7 O7+0.26 
' - z ' — 26 


358.5/363 




> 1.0 


2.19 


7 9n +0.30 
' -^ u — 30 


349.5/350 


MOS2 


> 0.3 


±.UU_q 32 


6.94i°; 2 ° 


381.9/401 




> 0.3 


2.19 


6.33lg;g 


407.7/402 




> 0.6 


2.19 


6.54l°; 24 4 


374.7/382 




> 0.8 


2.19 


a fi y+0.25 


354.6/369 




> 1.0 


2.19 


6 67 +0 - 29 

°'°'-n 29 


344.1/356 


pn 


> 0.3 


64+ ' 28 


n 77+O.33 
°-' '-0.33 


836.6/811 




> 0.3 


2.19 


°- ' D -0.19 


906.2/812 




> 0.6 


2.19 


fi 14+ - 28 

°- lq: -0.30 


818.8/754 




> 0.8 


2.19 


fi 49+ 031 

°-^ y -0.30 


729.0/714 




> 1.0 


2.19 


fi s^+o- 36 

O.80_o.35 


675.2/673 




> 1.2 


2.19 


7 01 + - 41 


636.6/633 




> 1.5 


2.19 


7 qo+0.65 
'•28_ 047 


582.9/572 



of this paper. To minimise these effects, we adopted the 
low-energy cutoffs derived above for the spatially-resolved 
analysis discussed below. 

5.2. Radial temperature profile 

We produced a radial temperature profile by excluding 
sources and extracting spectra in annuli centred on the 
peak of the X-ray emission. All spectra were binned to 3ct 
above background level (except the final annulus, which 
was binned to 2er) to allow the use of Gaussian statistics. 
We show the fitted spectra for annuli 5 and 9 in Figure |[ 
These spectra were fitted using the absorbed MEKAL 
model described above; we fitted separately the spectra 
from each instrument as well as maki ng a simultaneous 
MOS- 



-pn fit, as detailed in Table |B.l| . All the tempera- 
tures are consistent within the respective errors. It is also 
evident from Figure ^| that the form of each profile is sim- 
ilar. As a further test, we fitted the annular spectra with 
an absorbed MEKAL model with the absorption left as a 
free parameter. This produced profiles with, again, exactly 
the same shape, giving us high confidence in the form of 
the profile which we have derived. 

5.3. Projection and PSF effects 
5.3.1. PSF correction 

The PSF of XMM-Newton is a potential cause of con- 
cern, especially in the inner regions, where the bin sizes 
are small. To assess the effect of the PSF, we first calcu- 
late a redistribution matrix, F(i,j), where F(i,j) is the 
fractional flux in annulus i coming from annulus j. These 
redistribution factors were derived from our best model 
of the gas density profile, converted to emission measure 
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Fig. 4. XMM-Newton spectra of the cluster from an- 
nuli 5 (1.35' < 6 < 1.89', top panel) and 9 (5.13' < 
9 < 7.16', bottom panel). Black (red) [green] points: 
EPIC/MOSl(2)[pn] data. The EPIC spectra are back- 
ground subtracted and corrected for vignetting as de- 
scribed in Sec. || Solid lines: best fit isothermal model 
with parameters given in Table B.l. 



profile and convolved with the XMM-Newton PSF. The 
fractional contribution in each bin of the emission coming 
from the bin, as well as all inner and outer bins are plot- 
ted in Fig. ^. The PSF mostly affects the central regions 
and, above 2', the contamination from adjacent bins is less 
than 25%. 

We have a total of 30 spectra (10 annuli x 3 cameras) 
to be fitted with a model consisting of 10 MEKAL models 
(corresponding to the 10 'true' temperatures) absorbed by 
a common (frozen) absorption. The normalisations of the 
MEKAL models, for annulus i, are linked by the factors 
F(i,j), such as to reflect the contribution of each annulus 
j due to the PSF. In practice we ignore any contribu- 
tions at the less than 1% level. Each MEKAL model has 
6 parameters, which, together with the absorption, makes 
7. If we fit the 30 spectra simultaneously, the model has 
(6 x 10 + 1) x 30 = 1830 parameters. XSPEC can only han- 
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Fig. 5. The projected temperature profile of A1413. The bold 
black profile is the total total MOS +pn fit. For comparison we 
show the MOS 1 (red) MOS 2 (blue) and pn (green) separate 
fits. The error bars are la. 
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Fig. 6. Redistribution of the flux due to the XMM-Newton 
PSF: fraction of flux in each radial bin due to the emission 
of the bin (filled circles), as well as all inner (open circles) 
and outer bins (squares). 



die 1000 parameters, (even if most of them are frozen), so 
we have to find a way to reduce their total number. One 
way is to group the spectra, but for this to work the spec- 
tra in each group need a common normalisation. 

We typically find differences of ~ 20% between the 
normalisations of MOS and pn annular spectra (MOS1 and 
MOS2 normalisations are the same to within ~ 5%). We 
fit the global spectrum discussed in Sect. |5.1| to find the 
overall difference in normalisation between MOS and pn. 
We then checked that the annular fit results were the same 
when the global difference in normalisation was applied. 
This being so, we multiplied the pn annular spectra by this 
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Fig. 7. The effect of deprojection and the PSF. The projected 
temperature profile of A14f 3 (squares) compared with the pro- 
file obtained after correction for the PSF (diamonds) and the 
deprojected profile (triangles). Errors are la. 



factor to bring their normalisations into line with those for 
the MOS cameras. 

We now grouped the MOS and pn spectra of each 
annulus, giving 10 groups of 3 spectra, leaving us with 
(6 x 10 + 1) x 10 = 610 parameters, enabling a simultane- 
ous fit. We froze the metallicity of each MEKAL model at 
the best-fit value found for each projected annulus. The 
free parameters in the fit are then the temperature and 
normalisation of each MEKAL model. The resulting PSF- 
corrected profile is shown in Figure f7|(lcr errors). The PSF 
corrected results are entirely consistent with the projected 
temperature profile, with systematic differences of about 
half the la errors in the first 3 bins and smaller beyond. 
This result is not surprising, since the profile is relatively 
flat. Consideration of the PSF has a much smaller effect 
on the temperature profile of A1413 than for the bright 
cooling flow clust er A18 35 at z = 0.25 (Markevitch |2002 , 
Majerowicz et al. |2002[) . In contrast to A1835, A1413 dis- 
plays neither an extremely steep rise in the central gas 
density, nor a sharp drop in the temperature towards the 
center. As a result the contamination of central bins is first 
smaller: for A1835, more than 1/3 of the observed bright- 
ness at any radiu s is d ue to PSF scatterin g at smaller 
radii (Markevitch [2002 , Majerowicz et al 



2002j ) while 



for A1413 this contamination is already less than 25% at 
2', and decreases beyond (Fig. ^). Secondly, the smaller 
temperature gradient towards the center means that the 
redistribution biases less the temperature determinations. 

5.3.2. Spectral deprojection 

Another possible source of error in the derived profile 
comes from projection effects. 

A deprojected temperature profile was produced by 
first simultaneously fitting the MOS and pn spectra of the 
outer annulus with a MEKAL model absorbed by a fixed 



galactic column density. The spectrum of the next annu- 
lus inward was then fitted with a two-temperature model 
with the parameters of one of the models fixed to the 
best-fitting values derived for the outer spectrum. The 
normalisation of each fixed model must account for the 
volume within the outer shell projected along the line of 
sight toward the next shell inward. Furthermore, as the 
gas density profile is not flat, the normalisation must also 
account for this effect. We model the gas density profile 
using the parameters from the best-fit KBB model de- 
scribed in Sect. The normalisations are then adjusted 
by the emission weighted volume factors. This process was 
continued inward, adding a MEKAL model for each an- 
nulus, with the parameters of the outer annulus models 
frozen to their previously determined best-fit values. The 
abundances of the two outer annuli were frozen to the 
global value, so for these fits the free parameter is the tem- 
perature. For all other annuli both the temperature and 
abundance were free parameters. The deprojected temper- 
ature profile is shown compared to the projected profile in 
Figure [?]. In practice we find very little difference between 
the projected and deprojected results. The jump in the 
temperature of the ninth annulus is somewhat an artifact 
of the fitting process. In this case the software tries to 
compensate for the contribution of the low temperature 
found in the tenth (and first-fitted) annulus by putting a 
higher temperature in the subsequent annular bin. The er- 
ror on the tenth temperature is large, the contribution of 
the outer emission in the ninth bin depends on the actual 
cluster extent and thus the deprojected ninth temperature 
is probably more uncertain than found in this simple pro- 
cedure. Note, however, that the projected and deprojected 
temperatures agree well within the lcr errors. 

In summary, we find that neither a consideration of the 
PSF or projection effects substantially changes the form 
of the temperature profile. The profiles obtained by taking 
into account these effects are consistent with the projected 
profile, within the lcr errors. For all subsequent analysis, 
we thus used the projected profile. 

5.4. Modelling the temperature profile 

We now consider the scaled temperature profile, t(x) = 
T{r) /Tx, where Tx is the average temperature and x is the 
scaled radius, normalised to r2oo- ?"20Q is estimated from 
the average temperature Tx and the r2oo _ T relation of 
EMN96 at the cluster redshift. Tx is estimated by fitting 
the global spectrum, extracted from the [0.5' — 9'] region, 
i.e outside the possible cooling flow region (see below). 
We found kT x = 6.49 ± 0.15 keV (lcr error). Note that 
the temperature profile is determined up to ~ 0.7r2oo or 

~ ^500- 

We then modelled this projected temperature profile 
with a polytropic model: 



t(x) = r 



(7-1) 



(10) 
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Fig. 8. The projected temperature profile of A1413 with the 
best-fit polytropic model obtained with the central bin ex- 
cluded (full line). The best fit for the entire radial range is 
shown as a dashed line. 

where «h is the gas density profile given by the KBB 
model (Eq|), 7 is the polytropic index and To is the nor- 
malised temperature at x — 0. 

We fitted the profile with all (3, r c and £ parameters 
fixed at the values best fitting the surface brightness pro- 
file (Table ||), so that the free parameters for the fit are 
To and 7. When the whole radial range is fitted, we find 
r = 1.11 ±0.03 and 7 = 1.03 ±0.01, with X 2 = 7.62/8. If 
we then exclude the inner point, we find a better fit with 
to = 1.22 ± 0.03 and 7 = 1.07 ± 0.01, and X 2 = 4.45/7. 
The polytropic fits to the temperature profile are shown 
in Figure ||. 

We note that the fit is considerably better if the central 
point is excluded from the fit. The resulting polytropic 
profile rises to a peak in the centre which is not seen in the 
annular temperature determinations, which may lead us 
to believe that there is a small cooling flow (CF) at work 
in the very central regions. This possibility is discussed 
further in Sect. [| 

We further note that the derived value for the 7- 
parameter is very close to isothermal and moreover, is not 
very sensitive to whether the central bin is included in the 
fit. There is a tantalising hint that the temperature profile 
may drop further in the very outer regions, but the errors 
on this last data point are large enough that it is easily 
compatible with the derived polytropic model. 

6. Properties of the central gas 

Our best-fit polytropic model is an excellent description of 
the observed temperature profile barring the inner point, 
which is significantly lower. We examined if this temper- 
ature drop could be due to a cooling flow. 

The cooling time, the enthalpy of the ICM divided by 
the energy loss due to X-rays, is calculated using: 



tc 



2.9 10 10 vrs A /^% ( - 3 " H 3 )" (11) 
VI keV VlO" 3 cm- 3 / v ' 



Table 4. Multi-temperature and CF fits to the inner annulus. 
The F-test is computed against the fit for a single temperature 
absorbed mekal model 



Parameters 



IT 



2T 



CF 



kTi (keV) 


6.4 


6.9 


7.9 


kT 2 (keV) 




0.61 




Z/Z Q 


0.33 


0.35 


0.35 


M ( M yr" 1 ) 






58.9 


1003.2/854 


982.2/852 


978.3/851 


-fprob 




> 99.99% 


> 99.99% 



from Sarazin (1986). Using the central density derived 
from the KBB model fit (Table §), we find i coo i - 2.4 x 10 9 
yr, or about one quarter of the age of the Universe at 
the cluster redshift. This suggests that a CF should ex- 
ist. Furthermore, the cooling radius, defined as the radius 
where the cooling time is equal to the age of the Universe, 
is r coo i ~ 0.6' meaning that any CF should reside in the 
central bin. This is consistent with the observed temper- 
ature drop. 

We thus fitted the spectrum of the inner bin with more 
complicated models: 

— The sum of two MEKAL models absorbed by a com- 
mon column density fixed at the galactic value. When 
fitting, the second temperature is limited to be less 
than or equal to the temperature of the main compo- 
nent, and the abundances of the two components are 
tied together. 

— The sum of a MEKAL and a cooling flow model, again 
with a fixed common absorption. Here the abundance 
of the CF is tied to that of the thermal spectrum, and 
the upper temperature for the CF is limited to be the 
temperature of the thermal gas. 

The results are shown in Table ^. Both the two temper- 
ature and the MEKAL +CFLOW models are better fits 
than the single temperature model at the > 99.99% level. 
In addition, the MEKAL+CFLOW model is a better de- 
scription of the data than the two-temperature model at 
the 93% level. 

The secondary temperatures and CF properties are not 
well constrained. More sophisticated modelling is needed, 
preferably including RGS data, which is beyond the scope 
of this paper. 

7. Abundance 

The combined MOSFpn spectral fits also allow us to cal- 
culate a radial abundance profile, which can be compared 
with the profiles derived for cluster samples observed with 



BeppoS AX (I rwin & Bregman [2001] (IB01); De Grandi & 
Molendi [200l| (DM01)). Fig. § shows the profile derived in 
the conservative radial range where we have information 
from the Fe Ka line. 

The average abundance is 0.27 ± 0.03, more typical of 
the value found by DM01 for non-CF clusters (0.25±0.01) 



G.W. Pratt & M. Arnaud: An XMM-Newton observation of A1413 



11 



Radius(arcmin) 



0.7 

0.6 
^ 0.5 

M 
N 

— 0.4 

QJ 
V 

I 0.3 

c 

< 0.2 

0. 1 
0.0 



200 



400 600 
Radius (kpc) 



800 



1000 



Fig. 9. The radial abundance profile of A1413, derived in the 
radial range where the spectra have clear detections of the Fe 
K line. Errors are la. 



than for CF clusters (0.34±0.01). There is a clear decrease 
of abundance with radius. However, the profile exhibits a 
two-step behaviour, rather than a steadily declining pro- 
file, as found by DM01 and IB01. It remains more or less 
constant (a mean of 0.34, excluding the central bin) up to 
~ 400 kpc, dropping to another plateau (with a mean of 
0.22) beyond. IB01 found that inside 0.075r 20 o their CF 
subsample had an average abundance of 0.46 ±0.06, while 
the non-CF subsample had an average of 0.33 ± 0.04. The 
outer region abundances were 0.30 ± 0.02 and 0.24 ± 0.03 
for the CF and non-CF subsamples, respectively. The in- 
ner two bins of our observation correspond roughly to 
0.075r2oo- The mean error-weighted value for this region is 
0.34 ± 0.08 (3<7 errors, for direct comparison with IB01). 
Outside this region, the abundance value is 0.28 ± 0.09. 
Thus when the errors are taken into account, we can- 
not distinguish between the CF and non-CF subsamples 
of IB01. However, taken at face value, these figures ap- 
pear to suggest that the abundance profile shape of A1413 
displays characteristics intermediate between the CF and 
non-CF subsamples. This may be because A1413 appears 
to host only a modest CF. It must also noted that the 
abundance profiles of A1795 and A2142, both clusters with 
strong CF signatures (Peres et al. 
relatively flat abundance profiles, 
DM01. 



1998), appear to have 



as shown in IB01 and 
It has been suggested that both of these clusters 
are undergoing (or have undergone) mergers (Oegerle & 

1995D , which have presumably 



1994; Oegerle et al. 



Hill 

not been sufficiently strong to disrupt the CF, but which 
have effectively mixed the metals and thus flattened the 
radial profile. A study of cluster abundance profiles as a 
function of the strength of the CF signature and dynam- 
ical state would help to better understand the origin of 
cluster abundance gradients. 



8. Total Mass profile 



8.1. Calculation of the mass profile 

The mass profile is calculated under the usual assump- 
tions of hydrostatic equilibrium and spherical symmetry. 
The integrated mass profile can be calculated from the gas 
density, n g , and temperature profiles: 



M(r) = - 



kTr 
G/inip 



din n g 
rflnr 



rflnT 



din; 



(12) 



where G and m p are the gravitational constant and proton 
mass and /i = 0.597. 

If the gas density profile is described by the KBB 
model (Eq. [5]) , then the mass profile is described by: 



< R c 



where 



M{r) 



M(r) = - 



kr 
G/inip 

kr 2 
G/inip 



3Anr (2 «- 1) T(r) 



r 2 « + 
3/SYT(r) 



2« 



dT 
dr 



dT 
dr 



(13) 



, £, r c and /3 are the parameters of the KBB 
model, Pi n being linked to them by Eq. ^. 

The mass profile is calculated using the Monte Carlo 
method of Neumann & Bohringer ( 1995 ). The gas density 
profile parameters are fixed to their best fit values. This 
method calculates random temperature profiles within 
the bounds of the observed profile. We calculated 10000 
such random temperature profiles, using a window size of 
150kpc and a smoothing parameter of 0.1 keV. The final 
output is the mean mass profile and the corresponding er- 
rors for each data point in the input temperature profile. 
The errors are calculated using the 90% errors in the tem- 
perature profile and the standard deviation of the mass 
at any given radius. The resulting errors are calculated to 
correspond to la errors in the mass profile. 

The errors on the mass profile due to the error on the 
density gradient, dlnn g /<ilnr, are then calculated. As the 
parameters r Cj j n , £, r c and (3 are correlated, this error can- 
not be deduced directly from the errors on these param- 
eters. We used a method similar to the one described in 
Elbaz et al. (1995). For each considered radius, the sur- 
face brightness profile is fitted, considering d In n g /d In r 
estimated at this radius as a free parameter, instead of (3. 
The la error on this parameter is then classically derived 
from the \ 2 variation, the other parameters (normalisa- 
tion, r c .i n , r c and £) being optimised. Finally the errors 
due to the density and temperature profiles (derived from 
the Monte-Carlo method) are added quadratically. 

The resulting mass profile, with la error bars, is plot- 
ted in Fig. 0. 



8.2. Factors influencing the mass profile 

The temperature profile of A1413 is well determined out 
to ~ 0.7r2oo, and shows a gradual decline which is well 
described by a polytrope of index 7 = 1.07 beyond the CF 
region. Assuming such a polytropic description, the mass 
profile can be calculated analytically from the best fit gas 
density KBB model. The polytropic mass profile lies well 



12 



G.W. Pratt & M. Arnaud: An XMM-Newton observation of A1413 



10 



15 



04 

V 



on 



5 10 14 
o 
H 



Polytropic KBB model 




Isothermal KBB model 


Tcotrn^fTYi ill l"\ innrlf^l 

ISUlllCIlllcll L> 11HJUC1 








s? 




¥7 




/?' 




// 




// 

X/ 

: /// 


- 


// / 




/ / / 












/ / / , 





100 



Radius R (kpc) 



1000 



Fig. 10. The mass profile of A1413 derived from XMM- 
Newton surface brightness and temperature profiles. Data 
points: mass derived from the HE equation, using the 
Monte-Carlo method with the best fit KBB model for the 
gas density (Eq. ||) and the observed temperature profile. 
Errors bars are la and take into account both errors on the 
temperature and gas density profiles. Full line: mass profile 
derived from the HE equation and the best fit polytropic 
model with the KBB model for the gas density. Dashed 
line: same assuming isothermality. Dotted line: mass pro- 
file derived from the HE equation and an isothermal /3- 
model, fitting the outer gas density profile. 



within the errors of the Monte Carlo profile, except for 
the central point due to the temperature drop observed 
in the center (full line in Fig. [h]). We also note that the 
derived mass at large radii (r > 1.3 Mpc) lies at the lower 
range of the Monte Carlo mass. This is due to the drop of 
temperature in the last radial bin, the best fit temperature 
being below the polytropic value. 

In the classic approach, the gas density is described 
with an isothermal /3-model, in which the temperature 
profile is assumed to be isothermal and the gas density 
distribution is parameterised by a /3-model. In Fig. [w] 
(dotted line) we show the mass profile obtained using this 
approach, with the /3-model best fitting the outer part of 
the cluster and the average cluster temperature outside 
the CF region, kTx- Not surprisingly, the mass is greatly 
underestimated in the centre (r < i? C ut), where the gas 
density profile is more concentrated (higher gradient) than 
the extrapolated /3-model. If we instead parameterise the 
gas density using the best fit KBB model, the mass dis- 
tribution towards the centre is recover ed (Fig. [Uj[ dashed 
line). Beyond i? cu t the mass profile is slightly steeper than 
that derived from the true temperature profile, as ex- 
pected from the observed 7 value, slightly larger than 1. 



This comparison shows that the temperature gradient has 
a small but systematic effect on the derived mass profile. 

The mass profile is remarkably well constrained: the 
la error is less than ±5% below 1.4 Mpc and rises to 
~ ±18% at 1.8 Mpc. The temperature logarithmic gradi- 
ent is much smaller than the density logarithmic gradient 
(7% for 7 = 1.07), except in the very outer part, where the 
temperature gradient is both larger and the constraints are 
poorer. As a consequence (see Eq. |l^), except in this outer 
region, the error on the mass profile is dominated by the 
error on the gas density gradient (in the range 0.5% — 3%) 
and on the average temperature (2.3%). For the same rea- 
son the mass profile is very robust versus possible sys- 
tematic erro rs on the temperature profile. We have shown 
in Sect, b/2 that spectral deprojection or PSF correction 
do not have a significant effect on the form of the profile. 
One might also ask what effect the ellipticity of the cluster 
might have on the derived radial quantities. We extracted 
spectra in elliptical annuli and compared the projected 
temperature profile with that produced using circular an- 
nuli. All temperatures agree within the respective errors, 
and so we conclude that the cluster ellipticity is also a 
minor source of error. 

8.3. Modelling of the mass profile 



Navarro, Frenk & White (|1997| , NFW) performed high res- 
olution N-body simulations which showed that the density 
profiles of dark matter halos have a universal shape, re- 
gardless of halo mass and values of cosmological parame- 
ters. The NFW profile is given by: 



p(r) 



p c {z)5 c 



(r/r s )(l + r/r s y 



(14) 



where p(r) is the mass density and p c (z) is the critical 
density at the observed redshift, which, for a matter dom- 
inated Q. = 1, A = Universe is: 



Pc(z) 



(15) 



The parameters of the model are r s , a scale length 
and S c , a characteristic dimensionless density dependent 
on the formation epoch of the dark matter halo. S c can be 
expressed in term of the equivalent concentration param- 
eter, c: 



8c = 



200 



3 [ln(l + c) - c/(l + c)] 1 



(16) 



The radius corresponding to a density contrast of 200 is 
^200 = cr s . 



The NFW density profile varies from pnfw oc 



at 



small radii to Pnfw oc r 3 at large radii. As we are fitting 
the mass profile M(r), we use the integrated mass of th e 
NFW profile for the fit (e.g Suto, Sasaki & Makino |l998| ) : 



M(r) = Anp c (z)8 c rlm(r/r s 

X 

m{x) = ln(l + x) — 



(17) 
(18) 
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Fig. 11. The mass profile of A1413 fitt ed wi th the NFW 
profile (dotted line) and Moore et al. ( 1999 ) profile (full 
line). Bottom panels: residual between the data and the 
model. The best fit concentration parameters are c = 5.4 
and c = 2.5 respectively (see Table ||). The radius corre- 
sponding to a density contrast of 500 is indicated by an 



arrow. 



More recent, higher resolution simulations by Moore et 
al. (1999, hereafter MQGSL) suggest a profile described 
by: 



p(r) 



where 



5 C = 



p c (z)S c 



( r / r „)3/2 l + (r/rs )3/ 2 



lOOc^ 



ln(l + c 3 / 2 ) 



(19) 



(20) 



This is essentially identical to the NFW profile at large 
radii but is steeper near the centre (p m oc r -15 ). Again, 

Table 5. NFW and Moore et al. (1999) fits to the mass 
profile of A1413. Errors are ler 



Parameter 


NFW model 


MQGSL model 


c 


5.4 ±0.2 


2.6 ±0.1 


r s (kpc) 


401 ± 17 


845 ± 43 


r 2 oo (kpc) 


2169 


2221 


M 20 o( M ) 


8.9 x 10 14 


9.5 x 10 14 


X> 


8.76/9 


6.44/9 



as we are fitting the mass profile, we use the inte grate d 
mass of the MQGSL profile, given by (Suto et al. 1998): 



M(r) = 4TTp c (z)5 c r^m(r/r s ) 

2 ' • 3/2) 



m{x) = — ln(l 



(21) 
(22) 



The derived parameters from each fit are given in 
Table p| and the best fit models are compared to the data in 
Fig. [llj! We find that the data are extremely well described 
by the MQGSL profile across the entire radial range. The 
NFW profile can also be used to describe the data, but 
shows a small divergence at small radii (x 2 = 4.54 for 
the first 3 points). The radius corresponding to a density 
contrast of 500, computed from the data, is indicated by 
an arrow. There is a slight hint that the measured mass, 
M 500 = 7.7±Q 2 8 x 10 14 Mq, is higher than the MQGSL 
and NFW models (6.8 x 10 14 M and 6.5 x 10 14 M Q , 
respectively) around that radius (see Sect. !L3 for further 
discussion) . 

9. Discussion 

9.1. The shape of the temperature profile 

It is instructive to compare the projected XMM-Newton 
temperature profile of A1413 with the composite profiles 
found for larger cluster samples. The most extensive sam- 
ples come from ASCA and BeppoSAX data; these are, in 
order of publication: Markevitch et al. (1998; MFSV98), 
White (2000; WOO) Irwin & Bregman (2000; IB00) and De 
Grandi & Molendi (2002; DM02). The MFSV98 ASCA- 
derived profile is sharply decreasing, such that for a typ- 
ical 7 keV cluster the temperature drop is characterised 
by a polytropic index of 1.2-1.3. WOO finds that 90% of 
the cluster profiles in his ASCA sample are consistent 
with isothermality at the 3er-level. The IB00 BeppoSAX- 
derived profile extends only out to ~ 0.3r2oo and is 
flat or even slightly increasing. In contrast, the overall 
DM02 profile, from a larger sample of BeppoSAX obser- 
vations, is characterised by an isothermal core extending 
to ~ 0.2r2oo- Their CF subsample exhibits a temperature 
drop of a factor of 1.7 between ~ 0.2r2oo and ~ 0.5r2oo- 
The non-CF clusters exhibit a sharper temperature drop 
in the outer regions. DM02 suggest that an incorrect treat- 
ment of the BeppoSAX strongback may explain the dis- 
crepancy between their result and that of IB00. 

We have derived the projected temperature profile of 
A1413 out to ~ 0.7r2oo or ~ rsoo, in much finer detail 
than is possible with either ASCA or BeppoSAX . All in- 
dications are that A1413 is a relaxed cluster. 

Our data are compared to the DM02 and MFSV98 
composite profiles in Fig. Although each individual 
data point is (marginally) consistent with the typical re- 
gion defined by these composite profiles, there is an obvi- 
ous systematic difference in shape. The A1413 profile does 
not decline sharply like the composite profile of MFSV98, 
or the profile of DM02 beyond ~ 0.2r 200 - In Sect. |5~4|, 



we 



show that the polytropic model gives an acceptable fit to 
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Fig. 12. The projected scaled temperature profile of A1413 
compared to the composite CF cluster profile as obtained by De 



Grandi & Molendi 2002 (dot-dash lines are joining their data 
points plus or minus the la errors). The composite profile of 
Markevitch et al. L998 is shown as the shaded region (enclosing 
the scatter in their best fit profiles). The solid line is our best 
fit polytropic model (excluding the CF region). 



the data. Excluding the central bin, the 7 value (1-07) im- 
plies an almost isothermal temperature profile, and is not 
compatible with that found by MFSV98. It is very sim- 
ilar to that found by DM02 for CF clusters, but DM02 
reject the polytropic model on the grounds of poorness of 
fit, which is not surprising given the decline of a factor of 
1.7 in the temperature of their composite profile between 
0.2 and 0.5r 2 oo- We do not see a similar decline, and so a 
polytropic model is a good fit to these data. On the other 
hand, their best fit broken line model is a poor fit to our 
data: we find a \ 2 = 22.1/6 for their CF best-fit, and the 
fit is worse for their non-CF relation (\ 2 = 27.7/7). 

It is obvious that, given the extra radial range afforded 
by these XMM-Newton data, a flat or increasing profile, 
such as that of WOO or IBOO extrapolated to high radii, 
does not describe the A1413 data either. We emphasize 
again that the temperature gradient is modest: the tem- 
perature decreases by ~ 15(20)% between 0.1 r 2 oo and 
0.3(0.5) ^200- Fitting the temperature profile up to 0.3r2oo 
(i.e., excluding the last three temperature bins, and the 
inner bin) with a polytrope allows us to compare our pro- 
file directly with that of IBOO. We find r = 1.19 ± 0.03, 
7 = 1.06 ± 0.01, consistent with the value derived for the 
full radial range. This gradient is in agreement with the 
level of isothermality found by WOO and IBOO in that ra- 
dial range taking into account their errors, as well as that 



found by Allen et al. (2001a) from Chandra data below 
^2500 ~ 0-3 ^200- Our observation is also consistent with 
other XMM-Newton observations of nearby clusters, e.g. 
the slightly decreasing XMM-Newton temperature profile 
of Coma (10 % at 0.2 r 2 oo, Arnaud et al. |2001a|) and 



and found to be flat within ±10% beyond the CF region 
(0.1 r 20 o, Arnaud et al. |2001b| ). 



9.2. Shape of the total mass profile 

In Sec. |8.3j we showed that the NFW form can describe 
the total mass profile of A1413. However a slightly better 
agreement in the center is obtained with a MQGSL profile, 
derived from higher resolution simulations. 

With Chandra, it is possible to examine the central 
regions in great detail, at the expense of information at 
large radii. At present, it is unclear whether the NFW 
or MQGSL profiles provide the better description of the 
mass profiles derive d from Chandra observations. Allen, 
Schmidt & Fabian ( 2001b ) investigate several forms for 
the mass profile of RXJ1347. 5-1145, finding that both the 
NFW and MQGSL provide an acceptable fit, although the 
NFW profile is favoured in terms of \ 2 . Perhaps the high- 
est resolution examination of a cluster mass profile is that 
of Hydra A by David et al. ( 2001 ), who find p oc r -13 be- 
tween 30 and 200 kpc, which is intermediate between the 
NFW and MQGSL forms. The addition of a mass point 
from H Q observations leads them to favour the NFW pro- 
file, although the result is still consistent with the MQGSL 
result. 

One sticking point is the value of the concentration 
parameter from the NFW fit by David et al. ( |200l[ ). They 
find c = 12, which is 3 times that expected for a clus- 
ter of the mass of Hydra A. Interestingly, a similar value 
of c was found by Arabadjis, Bautz & Garmire (2002) 



from a Chandra study of EMSS 1358 +6245 On the other 
hand, the c parameters of Allen et al. (2001a) are better in 
agreement with the theoretical predictions. At large radii 
the NFW and MQGSL profiles coincide and we can com- 
pare the c value derived from our NFW fit, c = 5.4±0.2 
for a M200 ~ lO 15 ' 1 ^ 1 M cluster, to numerical simula- 



tions. Teyssier (2002) derived c parameters in the range 
4.9 — 9.5 for 5 clusters in this very mass range. The average 
c parameter derived by Eke et al. (1998) for hot massive 
clusters is c ~ 6. It must be noted that a relatively large 
dispersion on this parameter is expected from numerical 
simulations, with a l cr A( logc) = 0.18(50%) at a given 
mass (Bullock et al. 2001 ). In conclusion, we emphasize 
the excellent agreement in shape between the mass pro- 
file derived for A1413 and the theoretical expectations, all 
the more remarkable in view of the very small statistical 
errors on the profile. 

9.3. Normalisation of the M-T relation 

We now examine the normalisation of the mass profile. We 
will classically define Ms, the mass within a given radius 
r$, inside which the mean mass density is S times the 
critical density, p c {z) at the cluster redshift. For clusters 
obeying HE and self-similarity, the mass M$(T, z), scales 
with the cluster temperature and redshift as: 



the temperature profile of A1795, measured up to 0.4 7-200 



(23) 
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Fig. 13. The scaled mass profile of A1413 (data points), ex- 
pressed as a function of density contrast, 8, compared with 
the simulations of Evrard et al. (1996) (thick line). Filled cir- 
cles: SCDM cosmology. Open circles: ACDM cosmology. The 
solid and dotted lines represent the scaled best-fit Moore et al. 
(1999) profile and the scaled profile derived from the best fit 
polytropic model, respectively (SCDM cosmology). 



where h(z) is the Hubble constant normalised to its local 
value and M\q{8) is the normalisation at density contrast 
6 (here T\q is arbitrarily expressed in unit of 10 keV). 
The above relation is remarkably well verified by adia- 
batic numerical simulations, down to S ~ 200, for both 
SCDM and ACDM cosmolog y, wit h M w (5) independent 
of cosmology (e.g Mathiesen 2001). The variation of the 
normalisation Mio(S) with 6 reflects the (universal) clus- 
ter internal structure and is the same for all clusters for 
a given density contrast (although some scatter is present 
in practice). 



In Figure 13 we show the scaled mass profile of A1413 
defined as: 



M(S) 

where 
6 

with 

Pc(z) -- 



(1 m \ —3/2 
— I 
lOkeV J 

3M(r) 



47rr 3 p c (z) 



h 2 (z) 



3Hl 
8irG 



(24) 



(25) 



and h(z) = (l + z) 3 / 2 for the SCDM cosmology considered 
here. This profile can be compared with the normalisations 
Mio(8) derived from numerical simulations, allowing us 
to check the normalisation of the M$-T relation at differ- 
ent density contrasts. The results of EMN96 are indicated 
with a thick line. We also plot the scaled profiles corre- 
sponding to our best-fit MQGSL model (thin line), and 



to our best fit polytropic model (dashed line). The den- 
sity contrast is computed self-consistently for each profile 
using Eq. 

Both the data points and the MQGSL and polytropic 
model curves run parallel to the EMN96 profile down to 
5 ~ 600. This simply reflects the excellent agreement in 
shape of the A1413 profile with numerical simulations, as 
outlined above. Hp we ver there is a very significant offset in 
normalisation: Af (2500) = 7.3 ±0.2 x 10 14 M , compared 
with Mi (2500) = 9.95 x 10 14 M from EMN96. In other 
words, the predicted normalisation of the M-T relation 
lies ~ 36% higher than the observed value, in excellent 
agreement with the Chandra finding of Allen et al. ( [2001a , 
- 40% at 5= 2500). 

Below 6 — 600 the observed profile levels off, so that 
the data points seem to converge towards the EMN96 
predictions at small 8. At 6 = 500 M(500) = 1.8±^J x 
10 15 Mq, only 20% lower than the EMN96 value of 
A/io(500) = 2.2 x 10 15 M Q and actually consistent with 
it, especially if we also take into account the dispersion ob- 
served in the simulations (15%). If this effect is real, this 
would point to a fundamental difference of form in the 
total mass profile at large radii, i.e the real cluster dark 
matter density profile drops less steeply than the canonical 
r -3 law. The observed discrepancy in the normalisation of 
the M-T relation at high density contrast would thus be 
due to a flaw in the numerical simulations for the dark 
matter component. 

However, it is more likely that this level off of the ob- 
served profile is an artifact due to incomplete virialisa- 
tion. If there is residual kinetic energy due to infall, the 
HE equation applied to the observed temperature profile 
would over-estimate the true mass. Such an effect is ob- 
served in the simulations of EMN96, although its expected 
magnitude is somewhat smaller that we observe, ~ 10% at 
5 = 500 on average (but a large scatter exists). A further 
indication that this is the correct explanation comes from 
a comparison with the polytropic temperature model. Let 
us assume that the best fit MQGSL model indeed reflects 
the true total mass distribution. The fact that the mass 
profile derived from the best fit polytropic model closely 
follows this profile down to 5 = 200 (see Fig. |l3|), indicates 
that this model is a true representation of the thermody- 
namic state of the gas if it was in HE up to there. The 
drop in temperature in the last bin (5 > 550), as com- 
pared to this model, would thus be an direct indication of 
incomplete thermalisation. We also recall that there is a 
sudden drop in the surface brightness profile at 8 = 7.8', 
corresponding to a density contrast of 450 (computed with 
the MQGSL model) . This further supports our interpreta- 
tion: we might actually be seeing the expected drop of the 
X-ray brightness beyond the edge of the virialised (and 
hot) part of the cluster. Finally, incomplete equipartition 
between the electrons and the ions at the border of the 
cluster (Chieze, Alimi & Teyssier 1998) could also con- 



tribute to the low (electronic) temperature observed. 

The observed scaled mass profile depends on the cos- 
mological model via the function h(z), used in the seal- 
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ing, and the angular distance, cZa, used to convert angu- 
lar radius to physical radius. On the other hand the the- 
oretical normalisation A/ 10 (<5) appears to be insensitive 
to cosmology. We thus examined if a better agreement 
with the theoretical normalisation is obtained for the cur- 
rently most favoured ACDM model (f2 m = 0.3, Oa = 0.7). 
From Eq. [l^, Eq. [24| and JEq. |25|, the derived 5 value 
scales as (d,A.h(z))~^ and M as d\h(z), with h 2 (z) = 
O m (l + z) 3 + Q\. For a ACDM model (open circles in 
Fig. |l^), as compared to the SCDM model (filled circles), 
the data points are moved down and left along a line of 
slope 1/2 in the log-log plane, with 5 multiplied by 1.137 
and M multiplied by 0.938. The translation is modest and 
its slope is similar to the slope of the scaled mass profile 
around S = 1000 — 500, so the observed scaled mass pro- 
file remains essentially unchanged and the agreement with 
the theoretical curve is no better. 

In summary the good agreement between the mass pro- 
file shape and the numerical simulations, measured for the 
first time in the whole virialised part of the cluster, sug- 
gest that the modelling of the Dark Matter component 
is correct. However, the offset in the normalisation of the 
M-T relation suggests that some physics is lacking in the 
modeling of the gas. Several groups have studied the effect 
of non gravitational physics, like pre-heating or cooling, 
on the M-T relation (e.g. Loewenstein 2000; Bialek et al. 
1999|; Toz zi & Norman pOM] Babul et al. §00§ Thomas et 



al. |2002| ; Voit et al. |2002| j. A detailed comparison with 



numerical simulations, which would require a statistically 
representative sample, is beyond the scope of this paper. 
We simply note that pure pre-heating seems to have a 
small effect of the M -T relation in the high temperature 



range of A14 13 (L oewenstein 2000 Tozzi & Norman 2001 



Babul et al. P002[) an d that models including cooling (e.g. 
Thomas et al. 2002 ) seem to be more successful. We also 
emphasize that some care must be exercised when compar- 
ing theoretical predictions with these X-ray observations. 
The magnitude of the effect is not lar ge as compared to the 
typical difference of ~ 50% (Henry 200C) in the normal- 
isation derived by various groups using purely adiabatic 
simulations. The dispersion of the relation, observed for 
both simulated (~ 20%) and real clusters, require the use 
of statistically well controlled samples. Finally, a further 
ambiguity lies in the definition of the temperature. The 
temperature profile is not exactly isothermal (although 
our data suggest that the departure is small). Ideally we 
should compare the data from a given instrument using 
the spectral temperature estimated from simulation, af- 
ter full modelling of the plasma emission folded with the 
instrument response. The study of Mathiesen & Evrard 



(2001) indicates that the spectral temperature could be 
an underestimate of the mass-weighted temperature, by 
about 20%. Note that this effect would worsen the dis- 
crepancy observed above. All these systematic effects have 
to be well controlled, if we want to confirm the departure 
from the self-similar scaling and identify the physical pro- 
cess responsible for it. 
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Fig. 14. The gas mass fraction of A1413 as a function of over- 
density S. Data points: gas mass fraction obtained using the 
total mass derived from the Monte-Carlo method. Line: same 
using the best fit MQGSL mass profile. 



9.4. Gas distribution 

From our determination of the central gas density we have 
calculated a cooling radius of r coo i = 0'.6. Figure || shows 
that the gas density has begun to rise above the /3-model 
fit at about twice this radius. This suggests that the main 
driver of the central peak in the gas distribution is the 
cusp in the dark matter profile, not the CF. It is likely 
that the cusp in the dark matter profile acts as a focus for 
the gas and the CF. This has implications for CF mass 
deposition rates M deduced from excesses to the /3-model, 
in that any application to this cluster would result in a 
gross overestimate of M. 

The ratio of the X-ray gas mass to the total gravitat- 
ing mass is shown as a function of density contrast 5 in 
Figure |lj. The / gas rises slowly with increasing radius, up 
5 ~ 1000. After this point, the / gas either stabilises within 
the error if we use the total mass data points, or con- 
tinues increasing beyond if we used the best fit MQGSL 
total mass profile. A conservative estimate of the gas mass 

0.2 ± 0.02. This shows 



fraction at 5 ~ 500 is /, 
that the hot ICM is more extended than the dark matter 
distribution, as has been found in previous studies (e.g., 
David, Jones & Forman 1995 ), and is also seen in nu- 
merical simulations, and expected from purely dynamical 
reasons (Chieze, Teyssier & Alimi 1997 ). 

Assuming that the gas density profile follows the dark 
matter profile at large radii (within a factor of 2 be- 
tween r2oo/2 and 2r2oo) and a polytropic equation of 
state, Komatsu & Seljak (2001) derived an analytical so- 



lution of the gas in HE in an NFW profile. They obtained 
7 = 1.15 + 0.01(c — 0.5), or 7 = 1.14 for the concentra- 
tion parameter c = 5.4 derived for A1413, and a X-ray 
outer slope of /3 ~ 0.65 for a typical 6.5 keV cluster (their 
Fig. 14). The first assumption is roughly verified for our 
best fit profiles and not surprisingly, our derived values are 
in fair agreement with their model, although we obtain a 
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slightly but significantly lower value of 7 = 1.07 ± 0.01 
and a slightly higher /3 value. 

The value at 5 — 600, corresponding roughly to the 
virialised part of the cluster, / gas ~ 0.2, can be used to 
calculate the total mass density in the Universe, following 
the arguments that assume that the properties of clus- 
ters constitute a fair sample of those of the Universe as 
a whole (e.g., White et al. 1993| ). Assuming that the lu- 



minous baryonic mass in galaxies in A1413 is approxi- 
mately one-fifth of the X-ray ICM mass (e.g., White et al. 



1993|), and neglecting other possible sources of baryonic 
(f2t,/1.2/ gas ), where Qb is the mean 
0.0820 (O'Meara et al 



2001 



dark matter, Q m 
baryon density. For Qbh\ 

is the Hubble constant in units of 50 km s _1 Mpc~ x ) 
we obtain fi m = 0.3ih^Q 5 . 



10. Conclusions 

The main conclusions of this work may be summarised as 
follows: 



1 



of 



We have reported XMM-Newton observations 
A1413, a relaxed galaxy cluster at z=0.143. 

2. In a 2D (3 model analysis, we detect substructure to 
the south which does not appear to be interacting with 
the main cluster. 

3. Excluding the data from this region and all obvious 
point sources, we have measured the gas density and 
temperature profiles up to r 500 (corresponding to a 
density contrast S ~ 500, with respect to the criti- 
cal density at the redshift of the cluster). With the 
assumptions of HE and spherical symmetry, we have 
calculated the mass profile out to the same distance. 

4. The gas density profile is well described with a /3-model 
beyond ~ 250/i^q 1 kpc. We further parameterise the in- 
ner regions with a modified version of the /3-model (the 
KBB, Eq. H), which allows a more centrally peaked gas 
distribution. 

5. The temperature profile (excluding the inner point) 
is well described by a polytropic model with 7 = 
1.07±0.01. The decline is modest: a decrease of ~ 20% 
between 0.1r2oo and 0.5r2oo- 

6. The mass profile, derived from the HE equation, is de- 
termined with an accuracy of about ±5% up to r6oo 
and ~ ±18% at r5oo- It ca n be remarkably well de- 
scribed by a Moore et al. ( 1999 ) profile with a scale 
radius r s — 845 ± 43 kpc and concentration parameter 
c = 2.6 ± 0.1. An NFW profile also gives an accept- 
able fit but describes less well the central regions. The 
c values we find are in good agreement with those ex- 
pected from numerical simulations for a cluster of this 
mass. The Dark Matter modelling in these simulations 
is thus strongly supported by the excellent agreement 
between observed and simulated profiles. 

7. Beyond rgoo, the observed temperature and derived 
mass profiles begin to depart systematically from, re- 
spectively, the polytropic description and the Moore 



et al. (1999) profile. There is also a sudden drop of the 



surface brightness profile at r^Q. This suggests that 
the gas in these regions may not be in HE, and we may 
thus be seeing the outer edge of the virialized parts of 
this cluster. 

8. The offset in the normalisation of the Ms — T relation, 
with respect to the simulations of Evrard et al. (1996) 
is now confirmed to be ~ 40% across the entire ra- 
dial range up to r^oo (i.e., in the virialised part of the 
cluster). 

9. The gas distribution is peaked primarily as a result 
of the cusp in the dark matter profile. The gas mass 
fraction increases with increasing radius, to reach ~ 0.2 
at r 500 . 

We are now in a position directly to confront simula- 
tions with observations. The results are encouraging (the 
obvious validity of the modelling of the Dark Matter dis- 
tribution at large scale) but many questions remain. How 
peaked arc dark matter profiles? What is the relationship 
between between central dark matter cusps and CFs? Why 
are some studies finding unrealistic values of the concen- 
tration parameter? What is the source of the discrepancy 
in the M-T relation? 

The statistical errors on the observed quantities are 
now small enough so that we can determine in detail the 
intrinsic dispersion in cluster properties and systematic 
discrepancies with the classical self-similar model. To an- 
swer the above and other questions, a statistical sample 
of cluster properties would be of great help, preferably 
using Chandra to probe the central regions and XMM- 
Newton to determine properties at great distances from 
the cluster centre. Confrontation with numerical simula- 
tion is essential. The full range of observations, correla- 
tions between cluster properties, and detailed internal gas 
structure should be derived taking into account that they 
are viewed through a given instrument, so that we are able 
truly to compare like with like. 
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Appendix A: Cleaning of data for soft proton 
solar flares 

Various methods have been used to remove background 
flares from XMM-Newton and Chandra data sets. For 
XMM-Newton observations, perhaps the most simple is 
direct visual inspection of the binned high-energy [10-12] 
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Fig. A.l. Top: The MOS 1 [10-12 keV] light curve of A1413 
before (black) and after (red) cleaning for flares. Dashed lines 
show the ±3cr thresholds. Bottom: The Poisson fit to the his- 
togram of the light curve, from which the thresholds are cal- 
culated. The upper threshold is 14 cts/104s. 

keV light curve over the whole field of view, the adop- 
tion of a threshold level, and the exclusion of any inter- 
vals above the selected threshold (see e.g., Arnaud et al. 
2001). For Chandra data, an na clipping is used, where 
the average count rate in the [3 - 6] keV band over the 
field of view is calculated, thresholds are set depending on 
this value, and the light curve cleaned for any intervals 
where the thresholds are not met. 

It is known that the XMM-Newton quiescent back- 
ground level is variable by ~ ±10%, and so it is not possi- 
ble to set a rigid threshold level for flare rejection because 
of the risk of losing good data. The threshold level should 
ideally be dependent on the quiescent rate of the observa- 
tion in question. 

For the A1413 observation, we extracted the MOS [10- 
12] keV light curves in the field of view in 104s bins (chosen 
as this is an integral multiple of the frame readout time) . 
Similar light curves were made for the pn, but in the [12- 
14] keV band. We then made a histogram of each light 
curve and fitted this histogram with a Poisson distribution 



X s 



II 



(A.l) 



where the mean of the distribution, A, is the free param- 
eter of the fit. Following Poisson statistics, the error on 
the mean, a = \/X. We found that the Poisson distribu- 
tion was, without exception, a better fit than a Gaussian 
distribution. We then defined thresholds at the ±3tr level, 
and rejected any time intervals outside these thresholds. 
We show the Poisson fit and the original and cleaned 
light curves for A1413, which is a quiet observation, and 
for comparison, the observation of MKW9 (Neumann et 
al., in preparation), which shows several solar flares, in 



Figures A..1 and Figures A. 2 



Fig. A.2. Top: The MOS 1 [10-12 keV] light curve of MKW9 
before (black) and after (red) cleaning for flares. Dashed lines 
show the ±3<r thresholds. Note the y-axis is in log units. 
Bottom: The Poisson fit to the histogram of the light curve, 
from which the thresholds are calculated. The x-axis has 
been zoomed to emphasise the fit. The upper threshold is 19 
cts/104s. 

is not prone to the overestimation of the mean rate, the 
main problem with the na clipping method. 

Appendix B: Table 



This method is effective at finding the quiescent peri- 
ods, even for data strongly affected by flares. As such it 
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Table B.l. Radial temperature profile results. The spectrum of each annulus has been binned to 3a significance, except the final annulus, which has been binned to 2a. 
They are then fitted with a MEKAL model assuming a fixed absorption of Nh = 2.19 x 10 20 cm -2 . Errors are given at 90% confidence for one interesting parameter. Note 
that some abundance measurements have been frozen at the global value found for the cluster. 
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